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Abstract. 

In order to interpret observations influenced by dust and to perform detailed mod- 
eling of the observable characteristics of dust producing or dust containing objects, 
knowledge of the micro-physical properties of relevant dust species are needed. Labo- 
ratory measurements of cosmic dust analogous provides essential input for our under- 
standing of how dust particles can influences the dynamics and thermodynamics of the 
stellar atmosphere by their opacity. 

The formation of the dust grains influences the stellar atmosphere in two ways: In 
the gas phase chemistry, dust formation results in a depletion of certain elements, which 
influences the molecular composition of the gas, and consequently the corresponding 
opacities. On the other hand, dust grains have a rather high mass absorption coeffi- 
cient, which often may be comparable to the gas opacity or even exceed it. Due to its 
high opacity and the resulting radiative pressure, the dust has a strong influence on the 
structure of the atmosphere and the wind properties of AGB stars. 

Great care is needed when obtaining laboratory data as even a moderate variation 
of the different micro-physical dust values within the range expected for possible mate- 
rials has noticeable consequences for the interpretation of near-infrared colors of AGB 
stars. 



1. Galaxies care about AGB stars 



Cosmic dust plays a cr ucial role in the evolution of the Universe by assisting the 
formation of molecules (|Hirashita & Ferraral I2002T) . which are needed to cool down 
star f orming clouds for the very second generation of solar type stars to be able to 



form (Schneider et al. 2004). Furthermore dust is critical to the formation of plan- 
ets dJohansen et al. 20041) providing a clear link between understanding the where and 
when of dust formation and the epoch of planet formation on a cosmic scale. Finally, 
dust absorbs ultraviolet-optical light and subsequently re-emits it at infrared/millimeter 
wavelengt hs, strongly affectin g our view and hence our understandi ng of cosmic star- 
formation ( Steidel et all 120051) and quasars dUrry & Padovanilll995h . and our lack of 



understanding of dust remains the major uncertainty in observations of dark energy 
dDavis et al.ll2007h . 

The dust present in the Milky Way tod ay is mostly produced in the envelopes of 



evolv ed (age > 1 Gyr), low mass stars (e.g. jHofner & Andersen! 12007: Mat tsson et al. 



2008) . Observations of the presence of large amounts of dust in very early galaxies 
(Bertoldi et al. 2003|; iMaiolino et aDl2004h with ages less than 1 Gyr indicate that even 
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Figure 1 . The figure illustrates the minimum stellar mass (in units of M0 as a 
function of redshift for two different assumptions of when the onset of starformation 
occurred. Redshift 14 corresponds to an onset of starformation at 29 9 Myr after Big 
Bang when assuming a fiat universe with Hq = 70.2 and Qa = 0.725 dKomatsu et al.l 
l2010h . while redshift 12 corresponds to 371 Myr af Big Bang. 



though low mass stars seem to be the dominant dust producers in the present universe, 
it might not always have been the case. 

Theoretical studies indicate that under certain circumstances high mass stars which 
explode as supernovae may provide a fast and possibly efficient dust formatio n environ- 
ment in the early Universe (IBianchi & Schneiderll2007uNozawa et al.ll2007l) . However 
direct observational evidence for supernovae as a major source of dust is still missing, 
even in the loc al Universe, with the most reliable o bservations showing only a few 10~ 3 
M Q produced dRho et al.| [2008: Mei kle et al.ll2007t) . well below the cosmologically in- 
teresting limit of 0.1 M©. So currently there is no observational evidence for efficient 
enough dust production in relation to supernovae of any kind to explain the presence of 
dust in very early galaxies. Recent theoretical studies suggest that the high mass end of 
the asymptotic giant branch could be responsible for significant, rapid dust production, 
and calculated the approximate dust contribution that such stars could have, finding 
that they were possible viable sources for p art of the dust atz> 5 dValiante et al.ll2009[ : 



Dwek & Cherchnel l20Tfj| : lGall et aDl2010h . 
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2. The dust formation window in AGB stars 

The crude conditions for dust formation that needs to be met are relative low tempera- 



tures < 2000 K and reasonably high densities log n(H) > 10 s cm -3 (lSedlmayrj|1994l) . 
Beside this then the relevant condensation time scales also needs to be shorter than the 
time-period for which favorable conditions for grain formation hold during the evolu- 
tion of the considered object. 

In the interstellar medium the condensation time scale will typically be longer than 
the typical hydrodynamical timescale needed for dust grains to be able to condense out 
of the gas phase. It therefore seems that dust grains can only condense (nucleate) in the 
cool stellar extended winds of AGB stars, in the winds from giant stars and in supernova 
remnants a few years after the supernova explosion. 

The two main dust types that will form in the stellar environment will be carbon 
and silicates, although iron, corundum (AI2O3), magnetite (Fe204) as well as Si02 
grains are also likely to be present in some objects and to an extend where they need 
to be considered. However, the two main species will be carbon dust and silicate dust. 
These two dust types will typically not form in the exact same location due to the high 
bonding energy of carbon monoxide (CO). As CO is the most tightly bound molecule, 
then in an equilibrium situations all available carbon-atom and all available oxygen- 
atoms will form CO. 



2.1. Grain nucleation 

Nucleation is the first stage of the condensation process whereby a vapor transforms 
to a solid or liquid. This phase change requires some degree of supersaturation in 
order to drive the system through the relatively unstable reactive intermediates (clusters) 
between the atomic or molecular vapor and the macroscopic solid or liquid states. 

Presently only a few nucleation rates based on detailed calculations of relevance 
for AGB stars are available. One nice example is presented in these proceedings by 
Patzer et al. of nucleation studies of Ti C in the conditions of carbon-rich AGB star 
envelopes. For the M-stars Patzer (2007) have shown that it isn't necessarily the most 



abundant element which will form the nucleation precursors. 

Consequently, the nucleation rate which is a function of temperature, density and 
supersaturation (S) for a partic ular vapor is often calculated by either th e classical ho- 
mogeneous nucleation theory dBecker & Dorindll935l: Feder et ai1ll966h or by the re- 
lated scaled homogeneous nucleation theory (IHalelll986 ). 



The classical homogeneous nucleation theory was developed to describe the nucle- 
ation of volatile materials such as water, hydrocarbons or alcohols at relatively low lev- 
els of supersaturation (S ~ 1.1-5.0) and temperatures (~ 300 K). The theory describes 
the formation of critical nuclei in a supersaturated vapor by means of thermodynamic 
quantities. The essential basic assumption of this approach is that the properties of the 
clusters in the nucleation regime are given by the extrapolation of the bulk properties 
even into the domain of very small clusters or the interpolation of thermodynamic prop- 
erties between those of the molecules and the solid particles. With these assumptions 
both the thermodynamic functions such as entropy and enthalpy and the rate coeffi- 
cients describing cluster formation and destruction become simple analytical functions 
of the cluster size, which allow a straightforward calculation of the rate of formation of 
critical clusters. 
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A fundamental result of classical nucleation theory is the existence of a bottleneck 
for particle formation. The small unstable clusters which form at random from the gas 
phase have to grow beyond a certain critical size which corresponds to a maximum in 
the Gibbs free energy of formation and separates the domain of small unstable clusters 
from the large thermodynamically stable grains. The existence of such a critical cluster 
size also holds in more realistic theo ries of cluster for mation. However, a review of the 
available experimental literature by (Nut h et alj|20 00) shows that no experimental data 



exists to support the application of classical nucleation theory to the condensation of 
refractory vapors. Refractory vapors seem to condense out at different supersaturation 
ratio s than volati l e ma terials. 

ICherchnel fcOlCh describes an interesting way of dealing with the modeling of 
carbon- and oxygen-based grain nucleation using a chemical kinetic approach for car- 
bon and forsterite (Mg2 Si04). 

For carbon grains Keith & Lazzati (2010) have show that the nucleation rate is 



a function of the carbon concentration in the gas phase for a hydrogen-carbon gas at 
saturations 1.2 and 1.5 and temperature 4000 K. They find that for low saturations, 
higher H/C ratios drastically lower nucleation rates. 

In thin slices of presolar graphite grains extracted from carbonaceous chondrites 
it is seen that graphite grains contain refractory carbides (TiC, MoC as well as Fe-Ni 
metal). In some cases the sub-grains seem to be nucleatio n sites in other cases the y 
appear to have been captured by the growing graphite grain (Bern atowicz et al.1ll99ll) . 



2.2. Grain growth 

CO divides dust nucleation and growth into two different chemical paths 

• C/O > 1 Carbon chemistry (molecules: C 2 , CN, CH, C 2 H 2 , C 3 , HCN). 
Possible dust types: graphite (C), amorphous carbon (C), diamond (C), silicon 
carbide (SiC). 

• C/O < 1 Oxygen chemistry (molecules: OH, SiO, TiO, H 2 0, Ti0 2 , VO, ZrO, 
ScO, YO, LaO). 

Possible dust types: enstatite (MgSi03), olivine (MgFeSiOs), ferrosilite (FeSiOs), 
pyroxene (MgFeSi0 4 ), forsterite (Mg 2 Si0 4 ), fayalite (Fe 2 Si0 4 ). 

Carbon has the unique property that the atoms can form three different types of 
bonds through sp 1 , sp 2 (graphite) and sp 3 (diamond) hybridization. Amorphous carbon 
is a broad term covering materials which have a combination of the different bond types. 
Amorphous materials can show a whole range of different optical properties, related to 
the exact microphysical properties of the measured sample. Amorphous carbon is an 
illustrative e xample of this, as th e measured extinction can differ by a factor of 10 (see 
e.g. Fig. 2 in I Andersen et al.ll2003l) . depending on the detailed microphysical properties 
of the amorphous dust. 

Silicates are the most stable condensates formed from the abundant elements O, Si, 
Mg and Fe. Out of these four elements silicate grains form as silicatetrahedras (Si0 4 ) 
combined with Mg 2+ or Fe 2+ cations. In the crystalline lattice structures it is possible 
for the tetrahedras to share their oxygen atoms with other tetrahedras and thereby form 
many different types of silicates. The optical properties of these silicates all have reso- 
nances around 10-20 /mi, due to the Si-0 stretching and the O-Si-0 bending mode 
arising from the silicatetrahedras. Alignment of the tetrahedras may cause sharp peaked 
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resonances, whereas amorphous silicates will show a broad feature which can be seen 
as a blend of such sharp resonances. 

Example of oxide grain growth are e.g.: 

• Olivine: 

2xMg + 2(l-x)Fe + SiO + 3H 2 -> Mg^Fe^i-^SiO^s) + 3H 2 

• Pyroxene: 

xMg + (l-x)Fe + SiO + 2H 2 -> Mg. Y Fe ( i_. v) Si0 3 (s) + 2H 2 

It should be noted that most of the presolar silicates identified in meteorites as 
stellar condensates are not stoichiometric pyroxene or oliv ine and many of the iden- 



tified grains have sub-solar Fe/Si and Mg/Si c ompositions dNguyen et al.1 12010). It is 



also interesting to note that Hoppe et al. (2010) have found that presolar SiC grains (Z 



grains) from lower-than-solar-metallicity AGB stars are on average smaller than those 
from solar metallicity AGB stars. 



3. Dust spectroscopy - what is observed? 

A crystalline solid has a highly ordered lattice structure, with constant bond lengths 
and angles between atoms. In a solid, rotational motions are not possible and the 
vibrational-rotational transitions seen in gas-phase molecular spectra are replaced by a 
broad, continuous band at the vibrational frequencies. Due to the symmetry of the struc- 
ture, only a few of the possible lattice vibrational modes are optically active. Therefore, 
crystalline solids have only a few sharp features in the infrared. 
Observationally we find: 

• UV to IR absorption bands of interstellar dust (diffuse medium and molecular 
clouds). 

• Far IR (sub-mm) continuum emission from cold dust. 

• Photoluminescence bands of very small grains (extended red emission, aromatic 
IR bands). 

• IR emission (absorption) bands from warm (circumstellar) dust. 

To interpret the observed extinction related to a given type of dust material, it is 
necessary to determine the extinction efficiency (Q ex t)- The extinction is given as the 
sum of the absorption and the scattering (Q ex t = Qabs+Qsca) and is a function of the 
dimension-less size parameter x = 2a /A and a composition parameter, the complex 
refractive index of the material, m = n - ik. The problem is that of solving Maxwells 
equ ations with appropriate boundary condi tions at the grain surface. It was first done 



equations witn appropriate boundary conditions at tne gram surtace. it was tirst done 
by dLorenzl fl898l : iMiel [19081 : lDebydll909h and the method of solution is there often 



refereed to as Lorenz-Mie theory. 

Depending on the dust particles in question there are different ways of deriving a 
result. 

• For particles small compared to the wavelength the Rayleigh approximation for 
simple shapes or a distribution of ellipsoidal shapes (CDE) can be used. 
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• For particles large compared with the wavelength Geometrical Optics are needed. 

• For inhomogeneous particles effective media theory e.g. Maxwell-Garnett or 
Bruggeman which are "mixing rules" for evaluating the effective permittivity 
in terms of permittivites and volume fractions of each constituent material. 

• For nonspherical particles it is possible to use Purcell-Pennypacker method (DDA) 
or the T-matrix method. 

• BUT: Shapes are generally simple and there are problems for conducting materi- 
als or in strong absorption bands. 



4. What is measured in the lab? 



Several forms of laboratory data are needed to determine optical and spectroscopic 
properties of astronomical dust grains: Absorbance spectra, which can be converted to 
monochromatic mass absorption coefficients (/c) for direct comparison to observational 
spectra, under the assumption that the grain sizes and shapes of the laboratory samples 
are similar to those dispersed in space. Reflectivity spectra, which may be used to derive 
optical constants that are useful inputs for radiative transfer models and in analyzing the 
effect of grain morphology on spectral features. 



4.1. The optical constants n and k 

In order to evaluate the magnitude of size and shape effects, the complex optical con- 
stants of the material (or the dielectric function) as a function of wavelengths have to 
be known. The complex optical constants or dielectric constants are not constant! They 
are a function of wavelength and physical states such as physical density, crystallinity, 
grain size and orientation, etc. The integral Kramers-Kronig relations connect the real 
and imaginary parts of optical constants of a ma terial at a frequency poin t with their 



values over the whole frequency domain (see e.g. Boh ren & Huffmanlll983l) . 

Nearly all of the optical constants available have been measured either on bulk 
samples or on thin films. The reason for that is that the planar geometry of films or 
coatings is much easier to describe in deriving optical constants from a measured spec- 
trum than are irregularly shaped particles. However, even in the case of ideal bulk mea- 
surements, the determination of optical constants over a wide frequency or wavelength 
range is not a simple task. Since the material absorption in different spectral regions 
usually differs by many orders of magnitude. For the determination of k either trans- 
mission measurements on samples of very different thicknesses (from centimeter down 
to submicron scales) or transmission and reflection measurements have to be combined 
( Dorschner et al.|[l995l) . Many crystals show an anisotropy in their optical constants. 



In these cases, measurements with polarized light along the different axes of the crystal 
have to be carried out which require careful orientation o f the crystal and alignment of 
the polarizers. For more details see the recent reviews byEenning & Mutschkd (I2010T) 
and|Henning| d2010l) . 



4.2. Matrix effects 

The pellet technique is the classic laboratory measuring technique for determining ex- 
tinction properties of particle samples. With this technique a solid sample is mixed with 
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potassium bromide (KBr), cesium iodide (Csl), or polyethylene (PE) powder that have 
high transmission through certain IR wavelength ranges. The mixture is pressed at a 
10 Ton load which results in solid pellets of sizes 0.55 mm thick (1.2 mm for PE) and 
13 mm diameter which ar e easy to use for s pectroscopic analysis (e.g . iTamanai et al. 
20091 : IChihara et al1l2002l : IJaeser et alJ[l99l Koike & Hasegawal ll987h . 



The major advantages of the pellet technique are low cost, low sample consump- 
tion, longevity of the pellets in a desiccator and the exact amount of a measured sample 
is known which makes it easy to determine the monocromatic mass absorption coeffi- 
cient (a:). The main disadvantage is that there is the possibility of enviro nmental effects 



due t o the electromagnetic p olarization of the embed ding medium (e.g.lPapoular et al. 
19981 : iHenning & M utschkd l2000l : ISpeck et alJl2000l : IClement et al.ll2003l) . 



To compensate for the disadvantages of the pellet technique two other methods 
have been invented; th e diamond anvil cell ( e.g. iHofmeister & Bowevl 12006 ) and the 
aerosol technique (e.g. ITamanai et al1l2006l) . ITamanai et a l. (2006) demonstrated that 
the strong absorption peaks at approximately 9.8 and 11 /mi obtained by aerosol mea- 
surements for an olivinetype crystalline powder are shifted up to 0.24 //m compared 
with spectra obtained by the KBr pellet measurements but that weak s pectral features 



are af fected much less by the KBr medium effect. See e.g. Fig. 3 in ITamanai et all 
(2009) for an example of the shift between the aerosol measurement of Ti02 and the 
same sample embedded in KBr disclosing a peak shift from 13.53 //m to 15.61 /urn. 

Particles which have spherical or roundish shapes produce larger differences be- 
tween the spectra measured by the aerosol and the Csl pellet techniques as compared 
to the irregular shaped particles. 

While aerosol experiments and thin film of small particles in diamond anvil cell 
have the advantages of no matrix effects and provides a low-polarization environment, 
then they both suffer from that it is difficult to get a reliable estimate of the mass ab- 
sorption coefficient (k). 



4.3. Grain shape effects 

There is no general trends with grain shapes and grain morphology on the overall ex- 
tinction from grains. Based on knowledge of n and k it is possible to calculate analytic 
solutions of the light scattering problem for particles of arbitrary shape. However in 
many cases, spectra of irregular particles can be approximated by a suitably averaging 
over different ellipsoidal shape parameters. With these approximations it is possible to 
obtain simple expressions for an average extinction cross section. 

The spectroscopic measurements of dust particles in aerosol allow the investiga- 
tion of the actual morphology of the aerosol particles by filtering a nd subsequent scan- 



ning electron microscope (SEM) imaging. Tamanai et all ( 20061 120091) have studied 



the influence of morphological particle characteristics on dust spectra and find a strong 
dependence of the measured profiles of the infrared bands o n grain shape and agglomer- 
ation. Based on these measurements Mutschke et al. (2009) have come up with an inter- 
esting new method, where they investigate the ability of the statistical light-scattering 
model with a distribution of form factors (DFF) to reproduce measured infrared dust 
extinction spectra for particles that are small compared to the wavelength, i.e. in the 
size range of 1 /mi and smaller. 
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4.4. Temperature effects 

Almost all laboratory measurements are carried out at room temperature. There is 
therefore cause for concern about comparing predictions based on room temperature 
measurements with the cosmic environment, where the temperatures may be far lower 
or far higher. 



To investigate the possible temperature effects of carbonates iPosch et all (|2007) 
derived the infrared optical constants of calcite and dolomite from reflectance spectra, 
measured at 300, 200, 100, and 10 K, and calculated small-particle spectra for different 
grain shapes. They found in general that the grain shape variations dominate over 
possible temperature effects. Howev er, in the far-IR th e temperature effects seemed to 
become significant. However, when lKoike et all (|2006) investigated temperature effects 
on the peak position of the of the 49 and 69 /mi forsterite bands, the effect was smaller 
at shorter wavelengths. 

Niyogi et al. (these proceedings) presents the temperature and compositional ef- 
fects on spectral features of olivine minerals. All in all the conclusion is that importance 
of the temperature effect seem to differ for different minerals. 



4.5. /3-SiC vs. a-SiC 

It was previously suggested that the observed variability of the 11.3 /mi feature in 
some carbon star spectra could be attributed to different crystal structures of SiC (e.g. 
Speck et al.lfl997h . 



The crystal structure of SiC shows pronounced polytypism which means that there 
exist a number of possible crystal types differing in only one spatial direction. The basic 
units from which all polytypes are built are Si-C bilayers with a three-fold symmetry 
axis, in which the Si- and C-atoms are closely packed. Hence, each Si atom is tetrahe- 
drally surrounded by four C atoms and vice versa. 

If all shifts occur in the same direction, then an identical position of the bilayer 
in the projection along the hexagonal axis is reached after three stacking steps. The 
resulting structure is of cubic symmetry and because of the three-step stacking period 
this polytype is called 3C (C for cubic). Another name for this polytype, which is the 
only cubic one, is the often-used term yS-SiC. 

The other extreme is obtained, when the bilayers are shifted alternatingly in op- 
posite directions such that, in projection with the hexagonal axis, every other layer has 
the same position. The lattice is then of hexagonal type, and because of the two-step 
period the polytype is called 2H. The environment of each Si-C bilayer which has been 
produced by the alternating shifts is also called hexagonal. All other polytypes are built 
up by a characteristic sequence of cubic and hexagonal Si-C bilayers, for which the 
3C and 2H polytypes represent the limiting cases. For historical reasons any non-cubic 
polytype or mixture of polytypes is also designated by the term q-SiC. 

Detailed laboratory investigations byrMutschke et al.l(ll999h and lPitman et al.l (2008) 



show that there is no systematic dependence of the band profile on the crystal type for 
a- and /3-SiC. The observed variations seen in car bon star spectra is likely to be due to 
either different grain shapes dPapoular et al. 1998 ) a nd/or an increasing contribution by 



hydrogenated amorphous carbon ( Baron et al. 19871 : iGoebel et al1ll995l) . 
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4.6. Database of dust properties 

In order to evaluate the magnitude of size and shape effects, the complex optical con- 
stants of the material (or the dielectric function) as a function of wavelengths have 
to be known. A compilation of these quantities for astronomically relevant materials is 



present in the Heidelberg- Jena-St. Pe tersburg database at http://www.mpia-hd.mpg.de/HJPDOC 
(Hen ning et aDll999l : IJager et al.ll2003h . 



The database contains references to papers, data files and links to related Internet 
resources of measurements and calculations of optical constants for dust grains of as- 
trophysical relevance in a wavelength interval reaching from X-rays to radio. There is 
also information on amorphous and crystalline silicates, various ices, oxides, sulfides, 
carbides, carbonaceous species from amoiphous carbon to graphite and diamonds as 
well as other materials of astrophysical and terrestrial atmosphere interests. 

Acknowledgments. The author thanks Darach Watson for making the Fig. 1 Bar- 
bie colored diagram as well as for fruitful scientific discussions. The Dark Cosmology 
Centre is funded by the Danish National Research 
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